Mean-field theory of Bose-Fermi mixtures in optical lattices 
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We extend the results of [M. Lewenstein et al, Phys. Rev. Lett. 92, 050401 (2004)] and determine 
the phase diagram of a mixture of ultracold bosons and polarized fermions placed in an optical 
lattice using mean field theory. We obtain the analytic form of the phase boundaries separating the 
composite fermion phases that involve pairing of fermions with one or more bosons, or bosonic holes, 
from the bosonic superfluid coexisting with Fermi liquid. We compare the results with numerical 
simulations and discuss their validity and relevance for current experiments. We present a careful 
discussion of experimental requirements necessary to observe the composite fermions and investigate 
their properties. 



I. INTRODUCTION 

During the last few years, the intensive quest for the 
Bardeen-Cooper-Schrieffer (BCS) transition [J] in ultra- 
cold trapped Fermi gases has stimulated the interest in 
mixtures of fermions and bosons, from now on called 
Fermi-Bose (FB) mixtures. In these systems, sympa- 
thetic cooling techniques can be efficiently used to reach 
temperatures well down into the regime of quantum de- 
generacy 0- Recently, the rich physics of FB mixtures 
has become itself one of the central topics of the physics 
of ultracold gases. Various phenomena occurring in FB 
mixtures have been analyzed, including the phase separa- 
tion between bosons and fermions yjSJa], the existence 
of novel types of collective modes 0, 0, u3> Hi the ap- 
pearance of effective Fermi-Fermi interactions mediated 
by the bosons 0, HI E El j the collapse of the Fermi cloud 
in the presence of attractive interactions between bosons 
and fermions 113, or the effects characteristic for the ID 
FB mixtures flUll^. 

Recently, the possibility to load ultracold atomic gases 
in periodic potentials induced by laser standing waves 
(optical lattices) has attracted considerable attention, 
motivated partially by the resemblances of these systems 
to the solid-state ones 0|. The remarkable experimen- 
tal advances on atomic trapping and cooling on one side, 
and the possibility of manipulation of the interatomic 
interactions via Feshbach resonances 0] on the other, 
have opened the way towards a fascinating new physics, 
namely that of strongly correlated ultracold atomic gases 
in optical lattices. The most important result so far 
in this brand new research field has been provided by 
the achievement of the superfluid (SF) to Mott-insulator 
(MI) transition in bosonic gases [lj| Eg • 

The experimental success achieved in bosonic gases, 
has triggered the interest in the physics of BF mixtures 
in optical lattices, which under appropriate conditions 
can be well described by the Bose-Fermi Hubbard model 
(BFH) model 0, • A particularly interesting feature 
resulting from the BFH model is the possibility to pro- 
duce fermionic composites, formed by the pairing of a 
fermion and a boson (or a fermion and a bosonic hole) 
[l9l I20L l2lj L A Bose-Fermi lattice gas is somewhat sim- 



ilar to the Bose lattice gas of Refs. [15|,ll6j, although it 
exhibits a much more complex and richer behavior at low 
temperatures, as we have recently shown in Ref. [20|. 



In Ref. [2(j, we have discussed the limit of strong 
atom-atom interactions (strong coupling regime) at low 
temperatures, and have predicted the existence of novel 
quantum phases in lattices that involve the previously 
mentioned composite fermions, which for attractive (re- 
pulsive) interactions between fermions and bosons, are 
formed by fermion and bosons (bosonic holes). The 
resulting composite fermions present a remarkably rich 
phase diagram, and may form, depending on the system 
parameters, a normal Fermi liquid, a density wave, a su- 
perfluid liquid, or an insulator with fermionic domains. 



In this paper we significantly extend the analysis of 
Ref. [2(j. In Sec.[H]we present a description of the model, 
while in Sec. IIIII we remind the readers the main results 
of Ref. [2(|. In Sec. IIVI we determine the phase dia- 
gram of the system for arbitrary values of the chemical 
potential, the Fermi-Bose coupling, and the tunneling 
(hopping) amplitudes using mean-field theory j^. In 
this sense our results can be considered as a generaliza- 
tion of the seminal analysis of the Bose-Hubbard model 
by Fisher et al. |2^|. We obtain the analytic form of 
the phase boundaries separating the composite fermion 
phases from the phase consisting on a bosonic superfluid 
coexisting with a Fermi liquid. The method we use is a 
"traditional" mean field approach, which has a limited 
domain of applications, but has an advantage of being 
technically relatively simple. In the case considered in 
this paper, it can be rigorously confirmed using more 
sophisticated methods of functional integrals within the 
quantum field theoretical framework. In Sec. we com- 
pare analytical and numerical results, and discuss their 
validity and relevance for current experiments. Sec. lVIl is 
devoted to a careful discussion of the necessary require- 
ments to observe the composite fermions and investigate 
their properties in experiments. Finally, we present our 
conclusions in Sec. IVIII 
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II. BOSE-FERMI HUBBARD MODEL 

In this section we present a description of the model in 
question. Let us consider a sample of ultracold bosonic 
and (polarized) fermionic atoms trapped in an optical 
lattice, e.g. 7 Li- 6 Li or 87 Rb- 40 K. Due to the periodic- 
ity of the lattice potential, the single-atom states form 
energy bands, as in solid-state systems. If the tempera- 
ture is low enough and/or the lattice potential wells are 
sufficiently deep, the atoms can be assumed to occupy 
only the lowest energy band. Of course, for fermions this 
is only possible if their number is strictly smaller than 
the number of lattice sites (filling factor pp < 1). To 
describe the system under these conditions, we choose a 
particularly suitable set of single particle states in the 
lowest energy band, the so-called Wannier states, which 
are essentially localized at each lattice site. The system 
is then described by the tight-binding BFH model (for a 
derivation from a microscopic model, see Ref. [17|). which 
is a generalization of the fermionic Hubbard model, ex- 
tensively studied in condensed matter theory (c.f. [25|): 

Hbfh = - J2( jBb i b i + Jf/ ^ + h ' c -) 

(ij) 



E 



-Vrii{ni - 1) + Uniirii - p B ni - pftti., 



(1) 



where b\, bj, fj, fj are the bosonic and fermionic 
creation-annihilation operators, rii = b]b%, — fjfi, 
and hb and pf are the bosonic and fermionic chemical 
potentials, respectively. The BFH model describes: i) 
nearest neighbor boson (fermion) hopping, with an asso- 
ciated negative energy, —Jb (—Jf)', h) on-site repulsive 
boson-boson interactions with an associated energy V; 
iii) on-site boson-fermion interactions with an associated 
energy U, which is positive (negative) for repulsive (at- 
tractive) interactions |24j . 

For large hopping (3> U, V) the low temperature state 
of the system in 2D and 3D consists of a superfluid of 
bosons, which condenses in the zero momentum mode, 
while the fermions form a Fermi liquid. At exact zero 
temperature, the fermions would presumably form a su- 
perfluid state due to the p-wave pairing induced by the 
attractive boson-mediated fermion-fermion interactions 
0. In this letter we are interested in determining the 
boundary of this hopping-dominated fluid phase. This 
boundary is of primary experimental relevance, since it 
sets the possible conditions for the observation of the pre- 
dicted interaction-dominated quantum phases |20| |. In 
order to determine this boundary, we have to extend the 
analysis of Ref. |20| to the previously unexplored regime 
of Jb,f — U, V. To this aim we formulate the mean-field 
theory for bosons following the line of Refs. [22ll23| . We 
shall first consider the case of a small fermionic hopping 
(Jf — 0) and a fermionic filling factor pp not too close 
to 1/2. If pp — 1/2 the mean- field solution becomes un- 
reliable due to nesting effects and (in 2D) also due to van 
Hove singularities. 



III. COMPOSITE FERMIONS AND STRONGLY 
CORRELATED QUANTUM PHASES 



Let us first recall the results of Ref. |20j , which concern 
the limit of interactions much stronger than hopping. De- 
noting a = U/V, one obtains that in the case of vanishing 
hopping for p, — \p\ + s > a > jl — \fi] + s — 1, s holes (or 
alternatively —s bosons) form with a single fermion the 
corresponding composite fermion, which is annihilated by 

fi = y/{n - a)!/n!(6t)Vi ( vWfMf^/, )■ The 
inclusion of a small tunneling (^ U, V) as a perturbation 
leads to an effective Fermi Hubbard Hamiltonian: 



H, 



eff 



(ij) 



(i; , / J +l 1 .c.)+]f e// ^m 1 m J , (2) 

(ij) 



where m., = f} fi, —Jeff is the negative energy asso- 
ciated with the nearest neighbor hopping of composite 
fermions, and K e ff is the energy corresponding to the 
nearest neighbor composite fermion-fermion interactions, 
which may be repulsive (> 0) or attractive (< 0) 26]. 
This effective model is equivalent to that of spinless inter- 
acting fermions (c.f. |22ll2^| l. and, despite its simplicity, 
has a rich phase diagram. 

For K ef f > 0, and p F < 1 (or 1 - p F < 1), the 
ground state of H e ff corresponds to a Fermi liquid (a 
metal), and is well described in the Bloch representa- 
tion. In this limit, as discussed in Ref. [20|, the system is 
equivalent to a weakly-interacting (even for large K e ff) 
Fermi gas of spinless fermions (for pp <§C 1), or holes (for 
1 — Pf ^ !)• The weakly- interacting picture becomes in- 
adequate for pp — > 1/2, and for large A = K e ff/2J e ff, 
where the effects of the interactions between fermions 
become important, and one expects the appearance of 
localized phases. A physical insight on the properties of 
this regime can be obtained by using Gutzwiller Ansatz 
(GA) 20]. Defining A crit = (1 + ml) /{I - m% with 
m z = 2pp — 1, one can obtain that for A < A cr i t the 
ground state is a Fermi liquid, while for A > A cr i t it is 
a density wave with a period of two lattice sites. We ex- 
pect that the employed GA formalism predicts the phase 
boundary A crit accurately for p F close to 1/2. For the 
special case of half filling, pp = 1/2, the ground state 
is the so-called checkerboard state. For K e ff < 0, and 
when |A| <§; 1, and pp is close to zero (one), a very good 
approach to the ground state is given by a BCS ansatz, 
in which the composite fermions (holes) of opposite mo- 
mentum build p-wave Cooper pairs. The ground state 
becomes more complex for arbitrary pp, and for A ap- 
proaching — 1 from above. The system becomes strongly 
correlated, and the composite fermions in the SF phase 
may build not only pairs, but also triples, quadruples etc. 
The situation becomes, however, simpler when A < — 1, 
since the composite fermions group into domains, form- 
ing what one can call a "domain" insulator. 
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IV. "SIMPLE MAN'S" MEAN FIELD THEORY E = E + 2dJ B rip 2 + O(-0 4 ), where 



In order to observe the strongly correlated phases in- 
volving the composite fermions it is useful, if not nec- 
essary, to know where are their boundaries. Thus, this 
question has both fundamental and practical character. 
In this section we will answer this question using perhaps 
the simplest version of mean field theory. We stress, how- 
ever, that we have confirmed the results of this section 
using more sophisticated methods of functional integrals 
within the quantum field theoretical framework. 

In the following we shall consider the case of vanishing 
temperature, T — > 0, but still fc^T 3> Jf, in such a way 
that we can safely assume Jf — 0. In order to analyze 
the relevant intermediate case Jb — U, V , we employ a 
mean- field formalism [2^, HjJ . We shall analyze only 
homogeneous phases (it is possible to prove that within 
the mean-field limit supersolid phases are energetically 
unfavorable). We introduce the superfluid order param- 
eter ip = (bi) = (b\), for every site i. Then, neglect- 
ing higher order in the fluctuations, we can substitute 
b]h = i){b\+bi)-^ 2 . In this way, we can model the prop- 
erties of the complete Hamiltonian Hbfh by a sum of 
single-site Hamiltonians: Hmf = J2i(Hoi — i>Wi)> where 
H 0l = ^7ii(rii - 1) + Umrrn - p^Ftrii - + 2dJ B ip 2 , 

Wi = 2cIJb (bi + b\ J , and d is the spatial dimension. 
The mean- field Hamiltonian Hmf contains the same on- 
site terms as those of Hbhm and additional terms, which 
represent the influence of neighboring sites. 

As previously discussed, the ground-state of H oi con- 
sists of n(m) = n — sm bosons per site, where m is the 
number of fermions in the site, and s depends on the par- 
ticular region considered in the phase space. For a given 
fermionic filling factor pp , the ground-state presents at 
every site m — 1 (m = 0) fermions with probability pf 
(1 — pp)- Therefore the ground-state can be written as: 

I0o)(0o| = (1 — Pf)\ti = h, m = 0)(n = n, m = 0| 

+ pp \n = h — s, m = l){n = n — s, m = 1|. 

The zeroth-order energy is of the form Eq + 2d J sip 2 , 
where Eq = Eo(n, 0)(1 — pp) + Eo(n — s, 1)pf, with 

V 

Eo(n, m) — —n(n — 1) + Unm — fxgn — p,pm. (3) 

Due to the form of Wi , the lowest order correction intro- 
duced by the tunneling occurs at second-order: 

f |(n,0|>V|n,0)| 2 



1 + 2dJ B 



h + 1 



E (h,0) - E (n,0) 
\(n-s,l\W\n, 1)| 2 



Eo(n-s,l)-E (n,l) 



PF 



(4) 



e(n,0) e(n-l,0) 
h — s + 1 n — s 



e(n — s, 1) e(n — s ~ 1, 1) 



PF 



(5) 



where e(n, m) = ps — Vn + Um. If r > the system 
minimizes the energy by having ip = (normal phase), 
whereas if r < a nonzero ip (superfluid) is energeti- 
cally favorable. Therefore, the curve r = describes the 
boundaries between the phase with a superfluid bosonic 
gas, and the interaction-dominated phases. In Fig. ^ we 
have depicted the curves r = for different values of a, 
and different regions of the ps — Jb phase space. 



V. NUMERICAL RESULTS 

The analytical results of the previous section have been 
compared with numerical calculations based in the GA 
approach. Technical details concerning the application of 
GA for fermions in 2D and 3D are discussed in App. 
The GA method allows a generalization of the analytical 
results presented above to the case of finite Jf, and for in- 
homogeneous phases. In our numerical approach, we first 
evaluate the Bloch functions for the lowest energy band, 
and properly combine them to obtain the corresponding 
Wannier functions, largely localized at each lattice site. 
Using the Wannier functions it is easy to obtain the co- 
efficients Jb, J f , U and V [l^- We employ the GA 

wavefunction: \<p) = F]; Hn™o El,=o fn^rm K, m >, 
where n denotes the number of bosons and m the number 
of fermions, and neglect fermionic anticommutation rules 
between operators at different sites (for a detailed discus- 
sion on this point, see Add. lA"|l . We have checked that 
the value of the maximal bosonic occupation, N max , does 
not affect the calculations. The coefficients / n>m must 
satisfy Xmm \fn : m\ 2 = 1- In the following, we consider 
only homogeneous phases, although the calculation can 
be straightforwardly extended to inhomogeneous phases. 
Hence, the ground-state can be found by minimizing the 
energy on a single cell, while using periodic boundary 
conditions. The energy is minimized by changing the 
values of the coefficients f n ,m using a standard down- 
hill technique. The ground-state is calculated for a given 
number of bosons and fermions, and for the case {7 = 0, 
V — and a relative low lattice potential (large Jb). 

Starting from the chosen ground-state, we evolve in 
the phase space by adiabatically varying the parame- 
ters of the system. The time evolution of the f n .m 
coefficients is obtained by employing the proper mini- 
mization 



where W = J^VVj. Following the Landau argument 
for second-order phase transitions, one can easily write 



(ihd t — H) = 0. These equations deter- 
mine the dynamics of the /„.,„ coefficients in the phase 
space, and constitute the basis of what has been called 
dynamic GA (2^|. Since the total number of bosons 
and fermions is conserved during the time evolution, the 
chemical potential will be time dependent. This time 
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Figure 1: Phase diagram in 2D as a function of Jb/V, the 
fermion-boson interactions a — U/V, and the bosonic chem- 
ical potential jl = hb/V, for the case of Jf = 0. The 
lobes denote the analytical phase boundaries calculated us- 
ing our mean-field approach. For h — 1 < jl < n, and 
p, — n + s < a < p, — n + s + 1, the number of bosons, n, 
and the number of fermions, m, satisfy n + sm — h, and com- 
posites with one fermion and s bosonic holes (— s bosons) are 
formed. 
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Figure 2: Phase diagram as a function of the hopping J/V, 
and the bosonic chemical potential jl — hb/V, for a = 0.25 
and pf = 0.25. Thin solid lines: analytical results for Jp = 0; 
Bold solid lines: numerical results for Jf = 0; Dashed solid 
lines : numerical results for Jf = Jb- Phases A are formed by 
a Mott-Insulator phase for the bosons and a Fermi liquid for 
the fermions. Phases B are characterized by the formation of 
fermionic composites with one fermion and one bosonic hole. 

field approach provides only a qualitative picture. 



dependence can be evaluated by calculating two initial 
ground-states tpo(t = 0) -01 (t = 0) with, respectively, a 
number of bosons Nb and Nj, + 5N, where 5N -C Nb- 
We evolve the parallel trajectories while reducing Jb- 
The chemical potential can be then approximated as 

libit) « ((Mt)W(t)\Mt))-(Mt)\H(t)\Mt)))/SN. 

By launching various trajectories, we can explore those 
regions with an incommensurate total number of bosons 
plus fermions. Consequently, the trajectories do not en- 
ter into the regions of the phase diagram in which com- 
mensurate phases are expected. Therefore the expected 
lobular gaps are opened. As shown in Fig. [2 our numer- 
ical and analytical results are in good agreement. We 
must stress that the effects of the fermionic tunneling 
are not taken into account in the analytic calculations, 
but they can be easily included in our numerics. In 
Fig. |21 we also depict a case with finite fermionic tun- 
neling (Jb = Jf)- As expected, the lobes of the phases 
with composite fermions shrink due to the larger mobility 
of the fermions. 

Finally we would like to comment about the validity 
of the mean-field approach presented in this paper. In 
general, the mean-field approach is exact at dimension 
d = oo, and is expected to be reliable for n » 1, since 
the relative effects of fluctuations is then small. For the 
considered cases n ~ 1, in a general situation (not at the 
tips of the lobes) the upper critical dimension is d c = 2 
|23| . and therefore the mean- field approach is reliable. 
However, at the tips the transition belongs to a different 
universality class with d c = 3 p3]|. and hence our mean- 



VI. EXPERIMENTAL REQUIREMENTS 

Let us finally consider in more detail the necessary 
requirements to observe the composite fermions and 
to study their properties. Let us use the notation of 
Ref . [2fJ , according to which Roman numbers denote how 
many particles form a composite, and a bar over the num- 
ber indicates the composites fermion-bosonic hole(s). In 
the following discussion we shall assume Jp — Jb = J ■ 
We shall consider the regions of the phase diagram which 
are accessible at easiest: 

• Region I : In this region boson-fermion interactions 
are weak, and we deal with bare fermions. The 
hopping rate is of order J, i.e. relatively big in 
comparison to the effective interaction energy K e f f 
which is of order J 2 /V. 

• Region II : The boson-fermion interactions are at- 
tractive, and we deal with composites consisting of 
one fermion and one boson. Both the hopping rate 
J e f f , and the effective interaction energy K e f / are 
of order J 2 /V. 

• Region II : In this region the boson-fermion inter- 
actions are repulsive, and the composites consist of 
one fermion and one bosonic hole. As in Region II, 
the hopping rate J e //, and the effective interaction 
energy K e ff are of order J 2 /V. 

Note, that in the regions III, IV, etc. the effective 
interactions K e ff are of order J 2 /V, but the effective 
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hopping becomes of order J 3 /V 2 , J 4 /V 3 , and so on. The 
larger the composites are obviously the less mobile, and 
for this reason the observation of their physics is more 
difficult than for regions I, II or II. 

If we limit ourselves to the Bose- Fermi mixtures in the 
regime J < V, U, there will be then the following tem- 
perature regimes: 

1. T > \U\,V - This is the high temperature case. 
The system is a non-degenerated gas of fermions 
and bosons. 

2. T < \U\,V - The bosons enter the MI phase, 
and the fermions form composites (or remain bare 
fermions in the region I). Taking into account the 
fact that the Fermi energy for the effective parti- 
cles is provided by J e / / , we have to distinguish two 
cases: 

(a) When T > J, the composites constitute a 
non-degenerated Fermi gas in all considered 
regions I, II or II of the phase diagram. 

(b) For J 2 /V < T < J, in the region I a degen- 
erated gas is achieved, since T < Tp. Since 
T > K e ff, however, this gas will behave as 
a practically ideal gas, since the interaction 
effects will be masked by the thermal ones. 
Note, that in the other regions, T > Tp, be- 
cause Tp ~ Jeff - In another words, since the 
composites in the region II and II have a large 
effective mass, we shall deal in this regime 
of temperatures with a practically ideal non- 
degenerated gas of composites. 

(c) For T < J 2 /V, in all regions I, II and II 
a degenerated interacting gas of composite 
fermions is formed, since T < Tp and T < 

Keff- 

In order to discuss the feasibility of these different tem- 
perature regimes, let us introduce the following nota- 
tion: Let a = Vo/E rec be the optical lattice potential 
Vq in units of the photon recoil energy E rec , a z - the 
transversal oscillator length, and a - the atomic scatter- 
ing length (which, in principle, may be modified using 
Feshbach-resonance techniques). Approximating the on- 
site potential by a harmonic oscillator, we obtain that 
the interaction energy entering the Hubbard model is of 
the form 



V/E rec — 



According to the results of the present paper, one needs 
J/V < 0.05 in order to enter into the strongly correlated 
phases. That means, in turn, that in order to reach the 
temperatures in which a degenerated gas of composite 
Fermions may be observed, i.e. T < J 2 /V, one needs T < 

T . We recall that 
2.11 




/xK, whereas for Rubidium T rec = 0.17 /xK. Obviously, 
lighter atoms (Li) are more favorable. 

Of course, the other parameters may be adjusted in 
various ways. In particular, apart from Feshbach tech- 
nique, modifications of the transversal traps and lattice 
height can be used to control the atomic interactions. 
We have checked that, for instance, such manipulations 
in the case of 6 Li- 7 Li mixture, the value of a = U/V may 
be modified by a factor 4. This can be easily seen, since 

2* A* (1 + B) 



2.5 x 10~ 3 x y/a/2iT x (a/2a z )T rt 
for typical experiments, for Lithium atoms T ri 



(1 + VbTG)(1*B*D)' 

where A = clbf / a BB is the ratio of scattering lengths, 
B = tub I rap is the ratio of the masses, G = Vob/Vof is 
the ratio of optical lattice potentials, and D — uj z b/^ z f 
is the ratio of transverse trapping frequencies for bosons 
and fermions, respectively. 

For instance for 6 Li- 7 Li mixture B = 7/6, and taking 
G = 0.3, D = 0.3, we get a = 2o,bf/<ibb, we gain a 
factor of 2. Assuming some concrete values of the pa- 
rameters: clbb — 500ao (this value demands the employ 
of a Feshbach resonance), asF = 30ao (with ao being 
the Bohr radius), we get a = 3/25. If pp = 0.5 and 
oj z = 100kHz, one obtains then for the region I, U e ft = 
14.6nK, U BB = 51//K, U BF = 6.1/iK, J = 21.5nK. 

The above estimates imply that the previously intro- 
duced temperature regimes will be correspond to (1) 
T > 6 n K, (2) T < 6 fiK, (2b) T < 21nK (2c) T < 
15nK. From these estimations, and the current experi- 
mental state-of-the-art, we conclude that phase (2a), i.e. 
fermionic composites forming an ideal non-degenerated 
gas, should be relatively easy to achieve experimentally, 
since temperatures of the order of 1/iK are routinely ob- 
tained in current experiments. We would like to stress 
that such possibility alone could constitute a very inter- 
esting experiment by itself. The values to achieve degen- 
erated phases, of the order of lOnK, are more demanding 
experimentally, especially due to the difficulties to cool 
fermions imposed by Pauli-blocking. Nevertheless, this 
regime of temperatures is definitely not unrealistic, being 
at the border of what is nowadays achievable. The con- 
tinuous developments in sympathetic cooling of fermions 
and bosons allow to foresee that such a regime could be 
achieved routinely in optical lattices in the nearest future. 



VII. CONCLUSIONS 

In conclusion, we have completed the analysis of 
Ref. [2(| by including the tunneling effects within the 
mean-field theory. This analysis is of crucial experimen- 
tal importance since sets the conditions for the observ- 
ability of the interesting phases with composite fermions. 
We have also developed a numerical method which is in 
good agreement with our analytics, and additionally pro- 
vides information for those phases with finite fermionic 
hopping. Finally, we have analyzed the experimental re- 
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quirements to observe the predicted effects. In particu- 
lar, we have discussed different scenarios depending on 
the system temperature. 

We would like to stress that in this paper we have al- 
ways considered an homogeneous system (no additional 
external confinement). The analysis of finite temperature 
effects and inhomogeneities caused by the trap potential 
[30| is interesting in itself, and has been recently the sub- 
ject of investigations in Ref. 31] . Particularly interesting 
is that a random on-site chemical potential in the pres- 
ence of Fermi-Bose mixtures can lead to the achievement 
of various regimes of the physics of quantum fcrmionic 
disordered systems |32j . 
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Appendix A: VALIDITY OF SIMPLE MAN'S 
GUTZWILLER ANSATZ 

In this Appendix we discuss some technical details con- 
cerning the simple man's Gutzwiller ansatz for fermions, 
consisting in writing the variational wave function as a 
product of on-site Fermi operators, and neglecting the 
anti-commutation relations between Fermi operators in 
different sites. 

Formally, our variational approach is equivalent to re- 
placing the fermionic annihilation and creation operators 
fiifl by spin 1/2 operators cr(i), (i). This approach 



can be justified using the exact Jordan- Wigner transfor- 
mation in 2D and 3D ([H0II1, see also E3)- In thc 
case of 2D this transformation acquires the form 



S + U) = /jexp 



k^j 



S (j) = exp 



Szti) = fjfj - 5 



-i^2&rg(k,j)flf k 

k^j 

1 



(Al) 



fi, (A2) 



(A3) 



where arg(fc, j) is the angle between k— j and an arbitrary 
space direction on the lattice, which we choose as x. The 
fermionic Hamiltonian of Eq. (J2J becomes then 

Heff = J eff Y / S + (3y A ^ ) S-( i ) + h.c. 

+ K eff (ft® + (S z (j) + , (A4) 

where the "magnetic vector potential" 
A(i,j)= £[arg(M)-arg(fcJ)](S,(*0 + V2). (A5) 

k=^i,j 

The simple man's Gutzwiller approximation corresponds 
on this level to i) a variational ansatz for the ground 
state wave function in the form of product of on-site spin 
states, and ii) substitution of the "magnetic potential" 
by its average, which under assumption of mirror reflec- 
tion symmetry with respect to the lattice axes is zero. 
Similar approach may be applied in 3D, although in that 
case the 3D Jordan- Wigner transformation requires an 
extended Hilbcrt space and non-Abelian gauge transfor- 
mations l36l. 
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